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Abstract 

We propose an extension of the 4+1D SU(5) domain-wall braneworld of Davies, George and Volkas which includes the 
addition of a discrete A4 flavor symmetry. We show that lepton mixing and light Majorana neutrino masses can be generated 
from the additional A4 physics while at the same time sufficient parameter freedom can be maintained in the charged fermion 
sector to produce charged fermion masses and quark mixing naturally from the split fermion mechanism. Importantly, we show 
that the vacuum realignment problem typical of discrete flavor symmetry models of quark and lepton mixing can be solved 
by separating the appropriate flavons in the extra dimension, leading to exponentially sensitive suppression of the operators 
responsible for vacuum realignment. 

PACS numbers: 



I. INTRODUCTION 

The fermion mass hierarchy problem and the origins 
of quark and lepton mixing are three of the most out- 
standing problems of the Standard Model. The first con- 
cerns the 12 orders of magnitude spread amongst the 
Standard Model (SM) fermions, from the mass of the 
top quark down to the (upper bounds of the) masses of 
the neutrinos. The last two concern the Euler angles 
and CP phases of the unitary matrices which describe 
the relative strengths of the charged current interactions 
amongst the different mass eigenstates, which are the 
CKM matrix in the case of quarks and the PMNS ma- 
trix for leptons. For the CKM matrix, the Euler angles 
are approximately 12 = 13°, 13 = 0.2°, and 8 23 = 2.4° 
[I], while the CP phase is about S13 — 1.2. It is known 
that the PMNS matrix is approximately tribimaximal, 
with 12 = 34.06±Jif and 23 = 45° ±7.1° pQ. Further- 
more, recent experiments, in particular the results from 
Daya Bay and RENO [3] , support a small 613 for the 
PMNS mixing matrix while the nature of CP violation 
in the lepton sector is unknown. 

Extra-dimensional approaches to flavor are a com- 
monly used for the fermion mass hierarchy problem in 
particular. In extra-dimensional approaches to flavor, 
the hierarchy of Yukawa couplings to the Higgs is nor- 
mally generated from overlaps of fermion profiles which 
are exponentially dependant on the input parameters, 
an example of which is the split-fermion mechanism of 
Arkani-Hamed and Schmalz [3]. 

Given that there is no known reason for the apparent 
3 + 1-dimcnsional nature of our universe, it is interest- 
ing to consider whether there could be extra spatial di- 
mensions which are hidden with respect to the known 
four. In some extra-dimensional theories the hidden di- 
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mensions are compactified and miniscule. On the other 
hand, in braneworld scenarios, our universe is confined to 
a 3 + 1-dimcnsional brane in a higher dimensional space. 
Braneworld models first became popular with the advent 
of models with large extra dimensions which could solve 
the hierarchy problem, most notably the ADD model of 
Arkani-Hamed, Dimopolous and Dvali |S]. Later, Ran- 
dall and Sundrum showed that the hierarchy problem 
could be solved in the RSI model, in which there are two 
branes embedded in a slice of Anti-de Sitter space, with a 
UV brane on which gravity is localized and is strong and 
an IR brane on which the Standard Model (SM) fields are 
localized and gravity is weak [6]. The same two authors 
also proposed an alternative model, the RS2 model, in 
which the extra dimension of the warped spacetime was 
in fact infinite with gravity and the SM fields localized 
on a single fundamental brane [7]. For other foundational 
papers on extra-dimensional theories, see |B14T2~]. 

Domain walls were first proposed as a way to dynami- 
cally localize fields onto a lower dimensional subspace by 
Rubakov and Shaposhnikov [T3] • It has been shown that 
such topological defects can in fact dynamically local- 
ize gravity jl4j, yielding a dynamical realization of the 
RS2 model. This is attractive from an aesthetic point 
of view since the brane is no longer a fundamental ob- 
ject. In domain- wall braneworld models, chiral fermions 
are localized to the wall by Yukawa coupling to the sin- 
glet scalar field which contains the domain- wall kink, and 
scalars are similarly localized through quartic interac- 
tions with the kink. Gauge bosons are notoriously diffi- 
cult to localize, and the only known mechanism which is 
physically plausible is the Dvali-Shifman mechanism [15] , 
where the bulk is in the confining phase and respects a 
gauge group G which is broken by some scalar field on 
the wall to a subgroup H whose gauge bosons are then 
repelled by the bulk via confinement dynamics. 

Davies, George and Volkas proposed a phenomenolog- 
ically plausible domain-wall braneworld model in 4+1D 
spacetime in which the bulk gauge symmetry is SU(5) 
which is broken in the usual way to the Standard Model 
on the wall by an adjoint scalar field [16] . A consequence 



of the presence of this £[/(5)-breaking field in the back- 
ground along with the domain-wall kink is that the pro- 
files of the SM fermions are in general displaced from the 
center of the wall and are further split from their SU (5) 
multiplet partners, leading to a natural realization of the 
split fermion mechanism [4 . The authors then showed 
that in this particular model, the SM fermions could be 
split appropriately to generate the fermion mass hierar- 
chy naturally [T7]- It was shown in the same paper that 
generation of the quark mixing angles looked promising, 
and that proton decay could be suppressed by giving 
the colored Higgs scalar a large displacement from the 
brane. However, it was argued at the end of the paper 
that lepton mixing could not be naturally generated in 
the model while simultaneously generating the correct 
charged fermion mass spectra. This motivates us to use 
additional physics to solve the lepton mixing problem, 
and in particular to consider the addition of discrete fla- 
vor symmetries. 

Models with discrete flavor symmetries such as A4 rep- 
resent an interesting approach to explaining the quark 
and lepton mixing patterns. Models of this type were first 
explored in [T51E2] and references [5UJ HI] will prove to be 
relevant for our analysis. In the simplest model using A4 
in 3+ ID with just the Standard Model gauge group [21], 
typically the different mixing patterns are explained due 
to A4 being spontaneously broken to different subgroups 
in each sector: A4 —> Z3 in the charged fermion sector, 
and A<± — > Z2 in the neutrino sector. This is typically 
achieved by the addition of two A4 triplet Higgs fields 
which couple to different sectors and which attain dif- 
ferent vacuum expectation value (VEV) patterns. When 
this is done, the CKM matrix is found to be close to 
the identity and the PMNS matrix assumes a tribimaxi- 
mal form. When interactions between the two A4 triplet 
flavons are switched off, this arrangement is valid since 
the two non-aligned VEVs are both global minima of the 
potentials for each flavon. However, when interactions 
are switched on the two VEV patterns tend to align and 
thus the responsible cubic and quartic coupling constants 
have to be fine-tuned significantly to be small. This prob- 
lem is known as the vacuum realignment problem, and is 
typical of theories with discrete flavor symmetries which 
have extended Higgs sectors of this type. There are in 
general three ways to ensure that the troublesome inter- 
actions are suppressed. 

One is to make the theory supersymmetric so that the 
undesired terms are forbidden by holomorphy and renor- 
malization constraints on the superpotential [21] . An- 
other is to use additional discrete symmetries forbidding 
the interactions [52] . Yet another is to exploit the physics 
of extra dimensions, by localising the flavons on differ- 
ent branes or by splitting their extra-dimensional pro- 
files with very little overlap so that the interactions are 
naturally eliminated or very suppressed [23H25] . 

Given that extra-dimensional models have been very 
successful at explaining the hierarchy problem and can 
ameliorate one of the major problems of discrete flavor 



symmetry models, and that discrete flavor symmetries 
can reproduce realistic leptonic mixing patterns which 
can be difficult to produce in extra-dimensional models, 
the combination of the two approaches is quite attractive. 
There have already been many models in the literature 
uniting the two approaches, particularly with regards to 
the warped RSI scenario. Altarelli and Feruglio first pro- 
posed a model based on A4 with an SM gauge group and 
the flavons restricted to different branes [53] . There have 
also been models with GUTs [SUES! [27] as well as models 
with more complicated discrete flavor groups such as the 
double cover of A4, T' [55]. It was also shown by Kadosh 
and Pallante that the flavons could be put into the bulk 
to allow enough cross-talk between the flavons to gen- 
erate small quark mixing angles while at the same time 
maintaining the desired vacuum alignment |25j . It is also 
interesting to note that some 3+1D models with flavons 
in the 1' and 1" representations of A4, which our model 
also contains and were not previously considered in dis- 
crete flavor symmetry models, were proposed in [571 HH] • 
One of these [27] was also based on SU(5) x A4. 

In this paper, we extend the SU(5) 4+1D domain-wall 
braneworld model of Davies, George and Volkas [16] with 
the inclusion of a discrete A4 flavor symmetry group. 1 
We first set up the background configuration with a sin- 
glet scalar field forming a domain-wall kink and an ad- 
joint scalar field which attains a lump-like profile which 
breaks SU (5) to the Standard Model in the middle of the 
wall in order to facilitate the Dvali-Shifman mechanism. 
We then dynamically localize the required fermions and 
flavon Higgs fields embedded in appropriate SU(b) x A± 
representations to the wall via Yukawa and quartic cou- 
pling to the kink-lump solution respectively, and we give 
the forms of the profiles for the resultant localized SM 
components which are split according to their hyper- 
charges, yielding a natural realization of the split fermion 
mechanism. We show that the results in [TTJ with regard 
to the fermion mass hierarchy problem can be reproduced 
as well as quark mixing, neutrino mass squared differ- 
ences and a tribimaximal lepton mixing matrix from a set 
of 5D Yukawa parameters which are all of the same or- 
der of magnitude. In our model, it turns out the required 
scale of the breaking of A4 by the triplet flavons can be 
altered due to the fact that the Dirac masses for the neu- 
trinos can be suppressed by the split-fermion mechanism, 
and these scales can vary from the electroweak scale all 
the way up to the GUT scale. We finally show that split- 
ting the charged ^4-triplet flavon from the gauge singlet 
^-triplet flavon can exponentially suppress the interac- 
tions responsible for the vacuum realignment problem. 

In the next section we outline the basic background 
kink-lump configuration formed from a singlet scalar field 
which condenses to form the domain wall and an adjoint 
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scalar field which attains a non-zero vacuum expectation 
value on the wall breaking SU(5) to the Standard Model. 
Section Hill outlines both the fermionic and scalar matter 
content of our model as well as the SU(5) x A4 represen- 
tations to which they are assigned. Section |TV| and Sec. [V] 
then address the dynamical localization of the fermionic 
matter and the Higgs flavon scalars respectively. Section 
|VI| gives details of the electroweak Yukawa Lagrangian 
of the model and the forms of the fermion mass matrices 
that arise after the Aj-triplet flavons condense with the 
desired vacuum alignment. Our parameter fitting analy- 
sis yielding the desired fermion mass spectra, quark mix- 
ing, tribimaximal lepton mixing and the correct neutrino 
mass squared differences is given in Sec. |VII| In Sec. |VIII| 
we discuss our solution to the vacuum realignment prob- 
lem in our model, with the full flavon interaction poten- 



tials given in Appendix |A| Section IX is our conclusion 



II. THE BACKGROUND SU{5) DOMAIN- WALL 
BRANEWORLD CONFIGURATION 

In this section, we briefly cover the basic set-up of the 
model discussed in [TCI [TT] . In our RS2-like model, the 
extra-dimension y is infinite and the brane to which mat- 
ter is trapped is formed as a domain wall. To form a do- 
main wall, we need to introduce a singlet scalar field with 
a Z2-symmetric quartic potential with two discrete and 
disconnected vacua, and thus find a solution in which 
the scalar field interpolates between these vacua from 
y = —00 to y = +00. 

Gauge bosons are notoriously difficult to localize on do- 
main walls, and they cannot be treated in the same way 
as fermions and scalars. Instead we conjecture that the 
Dvali-Shifman mechanism |15j works in 4+1D. In gen- 
eral, the Dvali-Shifman conjecture states that if a gauge 
group G is confining in the bulk but is spontaneously 
broken to a subgroup H on the wall, then there should 
be a mass gap between the glueballs of G and those of H 
and thus localising iJ-bosons to the wall. For our model, 
G = 5(7(5) and H = SU(3) C x SU{2) I x U(1) Y , with the 
appropriate breaking achieved by an adjoint scalar field. 

Hence, the field content for the background is 



rj ~ 1, X ~ 24, 



(1) 



and the most general ^-symmetric potential for 
these fields may be written as 

V, 1X = (cry 2 - /4)Tr( X 2 ) + a V Tr( X 3 ) + \i[Tr{ X 2 )] 2 
+ X 2 Tr( X 4 ) + l( V 2 -v 2 ) 2 . 

(2) 

Finding the requisite background domain wall solution 
with the desired breaking of SU (5) on the wall involves 
finding a classical solution dependant solely on the extra- 
dimensional coordinate y of the Euler-Lagrange equa- 
tions. To find such a solution, we impose the boundary 
conditions 



where xi is the component of x correspond- 
ing to the appropriately normalised generator, 
dioff(2/3,2/3,2/3,-l,-l)V3/2\/5, proportional to 
hypercharge, and set all other components of x to zero. 
Numerical solutions exist for a finite region of parameter 
space. For the sake of simplicity and of yielding an 
analytic solution, we choose to impose the parameter 
conditions 

2/4(c - A) + (2cA - 4/A - c 2 )v 2 = 0, a = 0, (4) 

with A = Ai + ^jp. One can then show that the solution 
is 



r)(y) = v tanh(fcj/), xi(y) = Asech (ky), 



(5) 



where k 2 — cv 2 — jti 2 ,, A 2 = (2/j 2 , — cv 2 )/X. The solu- 
tion above is plotted in Fig. [T] Numerical solutions still 
exist for choices outside this parameter region, with the 
solution for r\ always being kink-like and that for xi be- 
ing lump-like. The stability of configurations such as Eq. 
[5] against, for example, formation of non-zero values for 
other components of x has been checked [T5] . 

One could also consider non-perturbative corrections 
to the background solution. Outside the wall, the physics 
of confining SU(b) is fundamentally non-perturbative 
and cannot easily be calculated. However, since the local- 
ization of matter occurs within a distance scale of order 
1/k we can ignore non-perturbative effects |16j . 




a 



extra dimension v 

FIG. 1: A plot of the profiles for 77 and xi arising in the 
background kink-lump solution. 

Before we proceed, it is worth mentioning that there 
are several scales in our theory. Firstly, since gauge and 
Yukawa interactions are non-renormalizable in 4+ ID, we 
must impose a cut-off energy scale Ajjy. In addition 
to the cut-off scale, we also have S'(7(5)-breaking scale 
A-su(5) — Xi(y — 0) ) the bulk SU (5) confinement scale A c 
and finally the inverse width of the domain wall Ajjw — 
k. As explained in further detail in |16) . the required 
hierarchy for the model to be viable is Ajjy > A-si/(5) > 
A-conf > Adw- 



v(y = ±°°) = ±«j xi(y = ±00) = 0, 



III. THE MATTER CONTENT AND A± REPRE- 
SENTATIONS 

We now need to introduce three generations of quarks 
(3) and leptons as well as Higgs fields embedded in represen- 
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tations of SU(5) x A4. As is usual for SU(5) grand uni- 
fied theories (GUTs), the lepton doublets and right-chiral 
down-type quarks are embedded into SU(5) quintets, 
while the quark doublets, right-chiral up-type quarks and 
right-chiral charged leptons are embedded into SU (5) de- 
cuplets. Right-chiral neutrinos are introduced as gauge 
singlets. We will not discuss the group theoretic proper- 
ties of A4 in this paper, see [TBU2T] for example. 

In addition to the representations under the gauge 
group and the discrete flavor symmetry we must also 
consider the transformation properties of the fields under 
the discrete Z2 reflection symmetry which ensures topo- 
logical stability of the domain wall. Since interactions 
which localize fermions to the domain wall are Yukawa 
interactions of the form rj^S^ and 77 has negative parity, 
we must have "i>^> — > — ^>^, which we can satisfy by 
choosing $ — > iT 5 ^ or W — » — iT 5l &. Scalars can have 
cither positive or negative parity. 

The representations of the fermions, denoted as (Ri, 
R 2 ), where R\ denotes SU(b) representation and R% de- 
notes the representation under A4, are chosen to be 



*5 



(5*, 1), *' 5 ~(5*, 1'), <~(5*, 1' 



(10, 1) for i -- 
N ~ (1,3) 



1,2,3 



(6) 



where N is an ^4 4 -triplet containing all three right-chiral 
neutrinos. Under the reflection symmetry, N — > — iT 5 N 
and all other fermions transform as "J — > iT 5l $>. 

For the Higgs sector, we require at least one Higgs 
quintet which contains an electroweak Higgs and a col- 
ored Higgs scalar, and some additional flavons as per 
usual in models with discrete flavor symmetries. Since 
the three fermion quintets ^' 5 , and tpg are in the 1, 
V , and 1" respectively, and since all the fermion decu- 
plets are singlets under A±, to form ^-invariant Yukawa 
interactions which generate charged lepton and down- 
type quark masses we similarly require a Higgs quintet 
under each of the A4 representations 1, 1', and 1". As all 
three generations of right-chiral neutrino are embedded 
into an A4 triplet, and since we must form Yukawa in- 
teractions involving this triplet and each of the fermion 
quintets to generate a Dirac neutrino mass matrix, we 
must have another Higgs quintet in the triplet represen- 
tation of A4. For the desired off-diagonal elements for the 
Majorana mass matrix for the neutrinos, we also need a 
gauge singlet Higgs scalar transforming as a triplet under 
A4. Thus, our field content for the Higgs sector can be 
summarized as 



*~(5*,1), 3'~(5M'), 3"~(5M") 
p~(5*,3), p~(l,3). 



(7) 



Under the Z2 reflection symmetry, all scalars except tp 
are chosen to have negative parity, while ip is chosen to 
have positive parity for reasons which will be discussed 
later in this paper. 



IV. LOCALIZATION OF CHIRAL FERMIONS 

To obtain a 3+1D effective field theory on the do- 
main wall and calculate the electroweak Yukawa coupling 
constants and masses arising in the effective theory, we 
must localize chiral fermion zero modes embedded in the 
fermionic fields described in the previous section to the 
background domain-wall configuration, since these zero 
modes will be our candidates as the chiral fermions of 
the SM. This means we must couple the fermion fields to 
the background fields r\ and \. 

Let's consider the right-chiral neutrinos first. Since the 
right-chiral neutrinos are embedded into A4 triplets and 
since they are gauge singlets, they couple to r\ only and 
the trapping interaction is simply 



Y vxN = -h ln {NN) ir] . 



(8) 



The 5D Dirac equation that results from this is thus 

iT M d M N + h lr , V (y)N = 0. (9) 

To examine the effective SM Yukawa interactions for 
the neutrinos in the effective 4D theory on the wall, we 
can ignore the Kaluza-Klein (KK) modes and consider 
only the localized zero mode of the field N. We thus 
can simply look for a solution of the form N(x, y) = 
/iv(y)^fl(^)i where /jv(y) is the zero mode profile and 
vr{x) is an ^-triplet of 4D massless right-chiral neutri- 
nos satisfying the ansatz 



i^d^ R (x) = 0, 
j 5 v R (x) = +v R (x). 



(10) 



Substituting this ansatz into Eq. [9j we find that the 
profile /jv(y) satifies the first order differential equation 



dfNjy) 
dy 



+ hir,v tanh (ky)f N (y) = 0, (11) 



which can be easily solved to yield 



; 1 n " 



/jv(y) = CVsech * (ky), 
= C N k^ sech' 11 " (y), 



(12) 



where y — ky, hi n — and CV = Cwk~^ are the 

non-dimcnsionalizcd extra-dimensional coordinate, back- 
ground Yukawa coupling constant and normalization fac- 
tor respectively. 

For the analysis in this paper it is convenient to al- 
ways work with dimcnsionless variables and functions, 
thus we define the non-dimensionalized profile fw(y) = 
fc~ 5 /jv(y) ancl normalize it to one to obtain the correct 
normalization for the effective 4D kinetic term for the 
zero mode i>r. Since any field increases in mass dimen- 
sion by half when the dimensionality of spacetime is in- 
creased by one, we non-dimensionalize the profiles for any 
effective 4D mode is the same way. 



4 



/„<>•) 



2.5 




2/0 


- \ 


11.5 


- \ 


/ L0 




/ °- 5 





files f&{y) 

~ + h*V tanh (ky) + h^A^- ~ sech (ky)] f* Y (y) = 0. 

From the above equation, we find that the non- 
dimensionalized profiles f^ Y (y) °f the left-chiral zero 
modes embedded in the quintets are given by 



bf Y (y) = ~K log [ cosh (y)j 



/if x Farctan [tanh(|)] 
(16) 



FIG. 2: The localized right-chiral neutrino triplet 
profile for the parameter choice hi v = 100. 



As can be seen in Fig. [ij the dimensionless profile fa is 
peaked about y — and decays exponentially away from 
the wall. Hence vr is indeed localized on the domain 
wall. 

Next, we consider the fermion quintets. We have one 
of each of the quintets in the 1, 1', and 1" representations 
of Ai, which means that due to A4-invariance, Yukawa 
interactions between different generations of the quintets 
and the background fields ij and x are forbidden. In 
this case, the coupling of each of these fermions to the 
background is given by 



Ay) 



(13) 



~ ~l • • — - T - 

for the field N, writing ^Py(x,y) = f£y( 
R = 1, 1', 1" and Y = +1, - 1, and ha 



Note the relative minus sign change between the interac- 
tions of rj with N and rj with the fermion quintets. This 
choice was made so that a positive , h' 5r . and h'^- cor- 
respond to the existence of left-chiral zero modes for the 
SM components of the respective quintets. 

To find the profiles of these left-chiral zero modes 
embedded in the quintets, we repeat the analysis done 

{y)ip^ Y {x) for 
having the zero 
modes ipf Y ( x ) satisfy the same ansatz as vr given in Eq. 
10 but with the second condition of right chirality re- 
placed with that of left chirality "Y 5 iP§y( x ) = ~' 1 1 ; 5y( x )- 
On substituting the ansatz into the 5D Dirac equation 
for 



iT M d M - h*r,(y) - yJl\h* x xi{y)]^ Y {x,y) = 0, 
we obtain the ordinary differential equation for the pro- 
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FIG. 3: The profiles for the localized lepton doublet L 
and right-chiral down- type quark Dr arising from an 
arbitrary fermion quintet tyf? for the parameter choice 

-100. 



h R 

n 5r] 



100 and hf x 



The consequence of the hypercharge dependency in the 
coupling to Xi is that the degeneracy of the profiles for 
the lepton doublets and the corresponding right-chiral 
down-type quarks is broken, and furthermore their pro- 
files are displaced from y — 0, meaning they are split. 
The coupling to r], h"L, roughly determines the widths of 
the profiles, while the higher the ratio h^ x /h§ v the more 
the localization centers of the profiles are displaced from 
y = 0. A plot of the profiles for a lepton doublet L and a 
right-chiral down- type quark Dr, in any representation 
R = 1, 1', 1" of A4, for the example parameter choice 



h R _ 
n 5ri — 



100. 



hf x = -100 is shown in Fig. [ij 



Finally, we consider the localization of matter embed- 
ded in the decuplets ^\ . Since all of the decuplets are 
in the trivial representation of A4, off-diagonal Yukawa 
couplings between ry (or x) an d different generations of 
are permitted, unlike the case for the fermion quin- 
tets. Therefore the most general coupling of the fermion 
decuplets to the background fields is 



(17) 
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The background Yukawa couplings /iio^ = (^ion) ancl 
^lOx = (^iox) can ^ e thought of as 3 x 3 matrices in the 
flavor space spanned by the initial 5D fields ^\ - -Due to 
the Hermiticity of the Lagrangian, both these matrices 
must be Hermitian and thus we can always choose a basis 
in which one of them is diagonal at the very least. To be 
able to pick a basis in which both hio^ and h\o x are di- 
agonal requires that they commute, ie. [/iio^, hio x ] = 0. 
In that case, since hio v — diag(h\ 0r .,hi 0r .,hi 0r .) and 
hio x = diag(h\ 0x , hl 0x , hf 0T] ), each generation of decu- 
plet "i>\ for i = 1, 2, 3 obeys the Dirac equation 



iT M d 



M 



3 Y 

Viv) - \lc^K 0x Xi(y) ¥ WY {x,y) = 0, 



forF 



5 2 
4 

3 



( 18 ) 

Writing ^[ 9Y (x,y) = P WY (yW WY (x) and a S ain requir- 
ing that "0ioF i s a left-chiral fermion which obeys the 
massless 4D Dirac equation, we easily find that the non- 
dimensionalized profiles fl 0Y (y) for the decuplet zero 
modes take the same form as those for the quintets, 



10Y 



WY 




log [ cosh (y)] 



(19) 



io x Y arctan [tanh(-)]. 



In the case that hio n and hio x do not commute, then 
we cannot find a basis in which all three flavors decouple 
and the flavor diagonal eigenbasis of the full operator 

^io>) + \J\\h\Q X is in fact y-dependant. The series of 
equations becomes a matrix differential equation which 
is very difficult to solve. This kind of scenario is called the 
twisted split fermion scenario and has been treated in the 
context of other models in [3TJ E2] . We do not consider 
this case in the analysis and for the sake of simplicity 
we assume that and h\o x commute in this paper. 
A plot of the profiles for the right-chiral electron-type 
lepton Eji, quark doublet Ql and right-chiral up- type 
quark Ur for the example parameter choice h\ 0n — 100, 

h\ 0x = 100 for some generation i is shown in Fig. |4j 

We have successfully shown that fermionic sector of the 
SM can be localized on the domain-wall brane. Now we 
must consider the localization of the Higgs scalars and 
show that electroweak symmetry breaking is possible. 



LOCALIZATION OF HIGGS FIELDS 




FIG. 4: Profiles for a right-chiral electron-type lepton 
Er, quark doublet Q and a right-chiral up- type quark 
Ur arising from an arbitrary fermion decuplet ^!\ for 
the parameter choice h\ 0ri = 100 and h\ 0x = 100. 



the flavons themselves, which do not contribute to the 
localization of the profiles at leading order. Hence it is 
sufficient to solely analyze the terms coupling the flavons 
to rj and x an( i t ne bulk masses of the flavons. 

For the quintet scalars $ = & R=1 7 $' = $ R=1 ', $" = 
$> R=1 and the ^-triplet p = <j> fl=3 , the localization 
potentials are easy to write down. They are 



+ 2A^ xl (^)t$«Tr( X 2 )+A^ x2 (^)t( x T ) 2 $^ 

+ Xg,R vx {^ R yx T ^ R V, for R=l, 1', 1", and 3. 

(20) 

The mode analysis for the quintets follows that for the 
Higgs quintet in the original SU(5) braneworld model 
described in [17]. Taking the ansatz 



Pry{v)4>r\{x), 



-rn R Y<i>RY(x), 



(21) 



and substituting it into the resultant 5D Klein- 
Gordon (KG) equation, one can show that the (non- 
dimensionalized) profiles for the modes of the Higgs quin- 
tets, p 7 R Y (y) satisfy a Schrddinger equation with a hyper- 
bolic Scarf potential, Vnsiy), which can be written as 



We now wish to localize the required Higgs scalars 
on the domain wall. This involves examining the Higgs 
scalar potential. As is typical of models with discrete 
flavor symmetries, we have an extended Higgs sector and 
the full Higgs potential is very complicated. However, 
most of these interactions are self-interactions amongst 



v H s(y)pRY(y) 



EryPHyIv), 



v HS (y) = a 2 



Pry 
dy 2 

\y + (B 2 ry - A 2 RY - A RY ) sech z (y) 
+Br Y (2Ar Y + 1) sech (y) tanh (y), 



a,.-* (22) 



G 



where Ary and Bry are denned as 



3Y 2 ~ 



Ary — ^ I ~ 1 + ( 2 [(^*«xi + -^- A *« X 2 _ A * R »7 ~ 



3Y 2 ~ 2 a - 3Y 2 



Y2 




fly 



3 y x 

5 2 A * n »?x 



2A 



fly 



1 



the bulk masses, KK mode masses and quartic coupling 
constants to n and y are non-dimensionalized as 



~2 



™fly 



r/v 



flV 



A$h X 2 — 



A; 2 ' 
X^r %2 A 2 
k 2 ' 



A<& H yl — 



X$,r vx — 



A: 2 ' 
X^r x1 A 2 
fc 2 ' 
A^h^iM 

_ k 2 ' 



(24) 



and Ery are the eigenvalues of the above potential, 
which in terms of the mode masses and fundamental con- 
stants in Eq. [24] are 

Ery = fn RY — p.%R — X$R n + A 2 RY - (25) 



It turns out the profiles of the n = modes of the quin- 
tet scalars fields have exactly the same form as for the 
chiral zero modes for the fermionic quintets described in 
the previous section, with the Ary playing a role anal- 
ogous to the hj^ and the Bry being analogous to the 

lie 



'5x> 



PRY(y) =Czn Y e- b **rW, 

b$n Y (y) = A RY log [ cosh (y)] 



2BryY arctan 



tanh(-) 



(28) 
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The eigenvalues of the hyperbolic Scarf potential are 
well known [33-35 ■ In the case that Ar Y > 0, it is 
known that there exists a set of discrete bound modes 
for n — 0, 1, L^4flyJ with eigenvalues 

E^ Y = 2nA RY -n 2 . (26) 

This gives the mass of the nth localized mode as 

ji%R + \q,R ri - (A B y - n) 2 . (27) 



~ 2 ~2 
m nRY 



The lowest energy modes which have the same SM 
charges as the electroweak Higgs doublet, the n = 0, 

Y = — 1 modes, are the ones we identify as our candi- 
dates for the fiavons of the effective 4D field theory on 
the wall. It should be noted that there are regions of 
parameter space where for a given 5D flavon field, more 
than one mode has a tachyonic mass. It is also possi- 
ble to choose parameters such that the modes for the 

Y = +2/3 components, which transform under SU(3) C , 
would attain tachyonic masses, which would be disas- 
trous since then SU(3) C would be broken on the wall. 
Thus, to maintain an unbroken SU(3) C while employ- 
ing electroweak symmetry breaking and for the sake of 
simplicity in the analysis of the electroweak sector, we 
choose parameters such that only the n — modes of the 
electroweak components of $ fl attain tachyonic masses 
on the wall while all modes of the colored Y = +2/3 
components attain positive squared masses. 



FIG. 5: The profiles of the localized electroweak Higgs 
<f) R and colored Higgs scalar <j) R for an arbitrary quintet 
scalar <fr R for parameters such that X^r^ — —7500, 



A$h x1 = 1500, 



A $R 2 = -75000, and \$r„ x = 2000 



From now on, we shall denote the profiles of the 
n = modes as PR=i,Y=-i(y) =Pw(y), PR=i>,Y=-i(y) = 
Pw(y), PR=i" ,Y=-i(y) = Pw"(y) and Jjjj=3,r=-i(s/) = 
Ppw(y) for the electroweak components and the same ex- 
cept with w replaced by c for the Y = +2/3 components. 
Similarly we will denote the corresponding 4D fields for 
these modes as (fr WjC (x), 4> w ',c'{x), <fi w ",c"( x ) and p w ,c(x) 
for Y = —1, +2/3 respectively. A plot of the profiles for 
the electroweak and colored scalar components of a Higgs 
quintet in any A4 representation is shown in Fig. [5[ 

Now we turn to the gauge singlet, ^-triplet scalar (p. 
The localization potential for (p is given by 



W v =/4(w)i + A ¥ 



1 ( W ) lV 2 + 2X vx ( W ) 1 Tr(x 2 ). (29) 



In a similar fashion to the analysis of the quintet 
scalars, writing down the corresponding Euler-Lagrange 
equation, then writing p{x,y) = J2P™(y)4>™{x) and 
□3+1^(0;) = — m 2 0™(x), we find that the modes of <p 
satisfy a Schrddinger equation with a well-known poten- 
tial, in this case the Poschl- Teller potential, 



dy 2 



d(d+l)tanh 2 (y)-d p%(y) = EX(y), (30) 
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where the parameter d and the eigenvalues E™ are given 
in terms of the fundamental constants and mode masses 
as 



(31) 



E% = rh 2 -fi 2 - A ¥ 



Given d > 0, there exists a tower of discrete local- 
ized modes with eigenvalues E™ = 2nd — n 2 , just as be- 
fore with the hyperbolic Scarf Potentials for the quintet 
scalars. Similarly, we identify the effective 4D gauge sin- 
glet, ^-triplet fiavon ip as the n = mode, and we 
choose parameters such that this mode is the only one 
which attains a tachyonic mass on the domain- wall brane. 
The mass squared for this fiavon localized to the wall is 
just 



and the profile for ipo, p Vo (y) is 



(32) 



(33) 



The lowest energy and tachyonic 4D mode of <p, tpo, is 
always localized at y = 0. A plot of the profile for the 
example parameter choice d = 500.00 is shown in Fig. [6] 
In contrast, the profile of the electroweak component of 
the other triplet, p, is in general not localized about 
y = 0. The natural splitting between the two A\ triplets 
will lead to solutions of the vacuum realignment problem 
since the splitting will naturally suppress the responsible 
scalar interactions. This will be covered more extensively 
in Sec. [yTTTl 



4r 



VI. THE ELECTROWEAK YUKAWA LA- 
GRANGIAN AND FERMION MASS TEXTURES 

Given we now have a set of localized modes for the 
fermions and flavons on the wall, we need to determine 
the electroweak Yukawa interactions in the model and 
the effective 4D mass textures after electroweak symme- 
try breaking. The full 5D electroweak Yukawa potential, 
Y E w, is given by 



Yew = hi(*s) c m ® + tii{%) c W w &' 

+ hpW^pN), + h' p W 5 {pN) v + h" p Wl{ P N) v , 
+ M(NN°)i + h v [(NN c ) 3s .^] x + h.c. 

(34) 

From this, we can deduce that the elements of the ef- 
fective 4D up-type quark mass matrix are given by 



(35) 



M% = ±h l l(<t> w ) / f Uk (y)f Ql (y) Pw (y) dy, 



4/iY (<t>w) / fu* R (y)fQi(y)Pw(y) 



where the electroweak Yukawa couplings are non- 
dimensionalized as til = h+k? . All electroweak Yukawas 
will be non-dimensionalized this way. 

For the down-type quark sector, the rows of the mass 
matrix are given by 



A I 



3] 



1 - 



V2 



fD R (y)f Q i{y)pw(y) dy, 
fD R (y)fQi(y)Pw"(y) dy, (36) 
fD R {y)fQi(y) P w'(y) dy, 



and similarly the columns of the electron-type lepton 
mass matrix is 



FIG. 6: The profile of the lowest energy mode, ipo, of 
the A4 singlet fiavon field (p for the parameter choice 
d = 500.00. 



Ail 1 - _ 
E V2 



Me 



h R {y)h{y)p w {y) dy, 
fE H (y)h'(y)Pw"(y) dy, (37) 
fEi(y)h"(y)pw(y) dy. 



The Dirac neutrino mass matrix comes from the in- 
teraction terms coupling the fcrmion quintet fields, the 
neutrino triplet and the Higgs (5*, 3)-field p. The form 
of this mass matrix is then dependent on the vacuum ex- 
pectation value pattern of the lowest energy mode of p, 
p w {x). For this paper, we require that p w takes the VEV 
pattern, 



(Pw) = {v p ,V p ,Vp). 



(38) 
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Under this alignment, the form of the Dirac neutrino 
mass matrix takes the form 



special case that m p 



m'' the left neutrino diag- 



D 



m p 


m p 


m p \ 


m' p 
m» p 


ujm' 

? n 
to m p 


uj 2 m' p 
uw!' p ) 



/ V5m p 







V3m' p 




(39) 



■U{u>), 



where the matrix U (uj) is given by 
where 




(40) 



h p v p 



m' p = h' p v p 



h" p v P 



fN(y)fL(y)p P w(y) dy, 
fN{y)fL'(y)P P w(y) dy, 
fN(y)fL"(y)p P w(y) dy. 



(41) 



The matrix U (w) is unitary, and such a factorization of a 
mass matrix where the mass eigenvalues can be arbitrary 
while the diagonalization angles are fixed is an example 
of form diagonalizability |36] . 

Lastly, the Majorana neutrino mass matrix is gener- 
ated from the bare Majorana mass term M(NN c )i and 
the Yukawa interaction term h v \[NN c )z s -p\ . The form 
of the Majorana mass matrix is obviously dependent on 
the VEV pattern of the lowest energy mode of ip, ip . 
Unlike for p w , we instead give the field po the VEV pat- 
tern, 

(¥> > = (0,tV,0). (42) 
This yields a Majorana neutrino mass matrix of the form, 



M, 



v : Maj 



where 



Clearly M 



M M v 
M 
,M V M 



fN(y)Pv>o(y) dy- 



(43) 



(44) 



v Maj i s diagonalized by the orthogonal matrix 



1 



'1 



-r 



v 2 l o 1 



(45) 



Noting that the effective left Majorana neutrino mass 
matrix is given by M L - M u ^ D M~^ Maj M^ Dl in the 



onalization matrix, V„l assumes the tribimaximal form 
V vL = U{tu)P, 

4= o \ 

(46) 



2 


i 





s/e, 


v/3 


UJ 






~V6 


^3 




J 2 


^ 2 




~Ve 


v/3 


V2 



The PMNS matrix which describes lepton mixing is 
defined by Vpmns — Vj* l V v l, where V e L is the left elec- 
tron diagonalization matrix. Hence, if we obtain mass 
textures for the charged leptons such that V e L ~ 1 we 
recover approximate tribimaximal lepton mixing which 
is favored by experiment. 

One also finds that in the special case m p — m' — 
m' p , the mass eigenvalues for the left Majorana neutrino 

mass eigenstates are M ^"^ , , and M ^2i • ^ or tnc 

purposes of this paper and for simplicity of analysis, we 
will choose parameters such that this condition is true. 
In the next section, we shall show that there exists a non- 
hierarchical parameter choice such that the Euler angles 
of the CKM matrix and the charged fermion masses are 
generated, while at the same time V e L ~ 1, yielding the 
correct lepton mixing patterns, and that the neutrino 
mass data can be satisfied in the case that m„ — m' 



VII. GENERATING THE FERMION MASS HI- 
ERACHY, THE CKM MATRIX AND LEPTON 
MIXING: AN EXAMPLE 

We now give an example parameter choice in which 
the fermion mass hierarchy, the Euler angles of the CKM 
matrix, lepton mixing and neutrino mass squared differ- 
ences are generated. With regard to the CKM matrix, at 
tree level the CP violating phase has to be put in by hand 
by giving the electroweak Yukawa coupling constants ap- 
propriate phases. We will for the sake of simplicity and 
clarity of the solution ignore this issue. Note that it is 
in principle possible for a CP phase in the CKM matrix 
to be generated from cross-talk between the quark and 
lepton sectors, since the PMNS matrix given in Eq. [46| 
contains non-zero Majorana CP phases. 

Before we begin the analysis of the fermion mass spec- 
tra, let us consider the constraints on the vacuum ex- 
pectation values of the gauge non-singlet Higgs fields 
coming from the masses of the W and Z bosons. As 
is typical with models with discrete flavor symmetries, 
our model contains an extended Higgs sector. The elec- 
troweak Higgs doublets which arise in the effective field 
theory on the wall are contained in the quintet Higgs 
flavons, <&, <£>', $" and p. To obtain the required W and 
Z boson masses, as is usual in a multiple Higgs doublet 
model, we must have that the sum of the squares of the 
vacua of these fields is equal to the square of the usual 
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SM Higgs vacuum expectation value of 174 GeV. In our 
specific model, this requirement is 



substituting the parameter choices 2 , 



^(4> w ) 2 + (4>' w ) 2 + (^) 2 + 



174 GeV, 



(47) 



where (<j) w ), (4>' w ) and (</>^) are the vacuum expectation 
values of the electroweak doublets arising from $, $' and 
<f>" respectively, with v p defined as before. For the sake of 
simplicity in our analysis, we will assume that the vacua 
obey (<t> w ) = {<P' W ) = (C> =v P = (174 GcV)/v^ = 
71.0 GeV. 



We first analyze the charged fermion sector. Since 
the charged fermion mass matrices are derived from 
Yukawa interactions with the Higgs fields and 
we first need to localize the electroweak components 
of these Higgs fields. For the purpose of this analysis, 
we choose parameters such that A R ^ Y =-i = 100 and 
B R . Y =-\ = —10, for R = 1, 1', 1". There exists a large 
region of parameter space spanned by the quartic cou- 
pling constants A^b^, A$r x1 , A$h x2 j A<j,H„ x an( l the 
Higgs squared masses n\ R which yields A R<Y =-i = 100 
and Br_y=-i = —10. Furthermore, a subset of this pa- 
rameter space is phenomenologically acceptable, namely 
that at least the lowest energy mode for the electroweak 
components attains a tachyonic mass on the wall, in- 
ducing electroweak symmetry breaking and for which all 
modes for the colored Higgs components have positive 
squared masses leaving SU(3) C intact. There also exist 
choices in this parameter region which satisfy all these 
constraints and which displace the colored Higgs well 
away from the wall, suppressing colored-Higgs-induced 
proton decay. An example from this parameter space is 



A^Hj) — —39725, Atf,j? x i 

10 , \$R r 
TO? 



-24396, A 



4.6886 x 



= 5189.8, and ^ R = 49700, which gives 
-o,b, y=-i = —25.0, m 2 lRY= _ 1 > for all n > 0, 
positive squared masses for all colored Higgs modes, and 
yields A R Y =2/3 — 9-86 and B RY=2 /3 = 64.6, which in 
turn produce a profile localized two characteristic lengths 
1/k to the left of the center of the domain wall for the 
lowest energy mode of the colored Higgs. 



Now that we have localized the electroweak Higgs 
flavons responsible for charged fermion masses after elec- 
troweak symmetry breaking, we can proceed to choose 
values for the background domain wall Yukawa coupling 
constants and electroweak Yukawa constants to generate 
the correct charged fermion mass spectra. We find that 





= 701.41, 


h 10 X 


= 304.55, 


h 2 


= 609.43, 


h 2 


= 263.73, 


h 3 


= 500.00, 


h 3 


= 188.63, 



(48) 



for the fermion decuplet background Yukawa coupling 
constants 





= 117.48, 


h 5x - 


-239.30, 




= 185.40, 


~K X = 


-274.64, 


K 


= 203.50, 


Kc = 


-254.82, 



(49) 



for the quintet background couplings, and 



hf 


hf 


hf\ 


/640.51 


504.28 


580.59 


hf 


hf 


hf) 


= 501.22 


481.66 


524.49 


hf 


hf 


hf) 


\129.87 


128.95 


431.03 



h 1 


h 2 


h 3 \ 


/800.00 


119.50 


119.00 


h'l 


h' 2 


h' 3 


= 119.50 


800.00 


119.50 


h'! 1 


h" 2 


h'l 3 J 


1 119.00 


119.50 


800.00 



(50) 



(51) 



for the charged fermion electroweak Yukawa couplings, 
we obtain the following mass matrices for the up-type 
quarks, down-type quarks, and electron-type leptons, all 
in units of MeV, respectively, 



Mr, 



/494.21 
1044.4 
\8037.2 



1044.4 
2389.3 
14474 



8037.6 
14475 
1.7139 x 10 5 , 



'4.6631 1.0987 3.2390\ 
M D = | 10.408 116.60 56.779 , 
,138.17 218.91 4245. 2J 



(52) 



(53) 



2 The appearance of dimcnsionlcss numbers of order 100 should 
not be taken to imply any loss of perturbativity in the Yukawa 
sector of the low-energy effective theory. These large numbers 
are produced by taking the ratio of a given dimensionful quantity 
with a smaller one, where the denominator has been chosen as a 
convenient reference point. To evaluate the validity of perturba- 
tion theory, one should look instead at the effective low-energy 
theory for the quantised field fluctuations about our classical do- 
main wall background solution. The low-lying modes are the 
discrete 3+ld KK excitations, and the strength of their Yukawa 
interactions are given by the effective 3+ld Yukawa coupling con- 
stants we have computed. By construction, the effective 3+ld 
electroweak Yukawa coupling constants have exactly the same 
values as those in the standard model, and are therefore pertur- 
bative. Since the underlying 4+ld theory is non-renormalisablc, 
the lack of UV completeness will eventually manifest through 
the appearance of strong coupling, and the likely loss of unitar- 
ity, above some high energy scale. This breakdown scale has not 
been computed, but it must be higher than the threshold for the 
continuum KK modes. The UV cutoff must be set low enough 
to avoid the likely pathologies, and some UV completion is, of 
course, implicitly assumed to exist. 
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Ah 



0.53478 
5.1121 x 1(T 2 
,5.6198 x 1(T 2 



15.847 
107.52 
13.762 



210.19> 

217.48 

1750.9, 



(54) 



The mass eigenvalues resulting from these mass matrices 
are 

m u = 2.49 MeV, m c = 1.27 GeV, m t = 173 GeV, 

(55) 

for the up, charm and top quarks, 

m d = 4.47 MeV. m s = 114 MeV, m b = 4.25 GeV, 

(56) 

for the down, strange and bottom quarks, and 

m e = 0.511 MeV, = 106 MeV, m T = 1.78 GeV, 

(57) 

for the electron, muon and tauon respectively. These 
results obviously are in agreement within error of current 
data on charged lepton masses [T] . 

With regards to quark mixing, the Euler angles, &x2i 
0j{j L , and ®% 3 L of the left up quark diagonalization ma- 
trix, Vjjl, derived from the up quark mass matrix in Eqn. 
[52] are approximately 

= 17.264°, 9 13 i = 2.6874°, 9 2 / 3 i = 4.8658°, 

(58) 

and the corresponding angles 6f^ L , Of^ L , and 9^ L of 
the down quark left diagonalization matrix, Vr> l , are ap- 
proximately 



Of/ = 4.2387°, 9^ = 1.8634°, 6^ = 2.9746 



DL 



(59) 



from which one derives that the Euler angles 8j" 2 J 



Q.CKM 



&\ 3 KM , and 9 23 of the CKM matrix, given V C km 



V UL V, 



DLi 



are 



r^CKM 
u 12 



= 13.0°, 8 



CKM 
13 



= 0.201°, e 23 



CKM _ 2 ggo 

( 60 ) 

These results are in agreement with the current data on 
quark mixing angles pQ. 

Finally, the Euler angles 9ff, 9ff, and Q%\ of the left 
electron diagonalization matrix, V e £, are indeed small 



6g = 7.32 x 10~ 2 °, 9^=4.15x10" 



9|f = 0.925°. 



(61) 



Thus, the left electron diagonalization matrix is indeed 
very close to the identity matrix, which means in choos- 
ing parameters such that m p = m' = m' p , we get a tribi- 
maximal PMNS matrix generated entirely from the left 
diagonalization matrix for the neutrinos. The only re- 
maining things to check with regard to the mass fitting 
are the existence of parameter choices for the right chiral 
neutrinos generating the desired mass spectra and if, in 
the case m p = m' p — m p , the experimentally measured 
neutrino squared mass differences can be satisfied. 

In the analysis of neutrino masses, we give two exam- 
ple parameter choices for the non-dimensionalized cou- 
pling of the right-handed neutrinos to the domain-wall, 



hi n , the bare Majorana mass M, the non-dimensionalized 
Yukawa coupling h v describing the strength of the inter- 
action between N and the gauge singlet Higgs field <p, as 
well as the quartic coupling parameters determining the 
localization of the lowest energy modes for the Aj-triplet 
flavon fields p and tp. The Dirac masses m p , m', ml' are 
sensitive to the splittings between the fields p, the lep- 
ton doublets and N, and since we still have the freedom 
to shift p, we can naturally control the scales of these 
masses: they can be made to be at the same order as the 
top quark mass or they can be made to be very light. On 
the other hand, since the right-handed neutrino triplet 
N and the flavon tp are always localized about y = 0, 
the overlap integral J /jv(z/)ffyo(2/) dy on which the mass 
scale M v is dependent is always naturally of order 1, so 
naturally we expect that M v is of the same order as the 
bare Majorana mass for N. 

In our first examples, we give parameters such that the 
neutrino Dirac masses are of the order of the electron 
mass. For the parameter choice 



h lv = 100.00, 
\ p = 533.27, h' p = 267.97, 



A 



pY=-l 
1 



B P Y= 



^«=3,y=-i 

-B»=3 V=—l 



>R=3,Y 

we get the Dirac masses to be 



h' p = 777.57, 
= 500.00, 
= 380.00, 



(62) 



m p = m! p = m" p = 0.100 MeV. 



(63) 



One can then show that to satisfy the neutrino mass 
squared differences Am 2 2 = 7.59 x 10 -5 eV 2 and Am 2 3 = 
2.43 x 10~ 3 eV 2 from the PDG [J, we must have 



M = 2.86 TeV, M w = 2.26 TeV, 



for which the neutrino masses turn out to be 



(64) 



mi 



m 2 



ma 



-3to 2 



M + M, r , 



5.86 x 10" 3 eV, 



-3m 2 



M 



-3m 2 



M - M, 



= 0.0105 eV, 



0.0504 eV, 



(65) 



Because the localized mode of tp, tp is a gauge singlet and 
A 4 -triplet and thus decoupled from the Standard Model, 
at least at tree level, we have the freedom to choose the 
scale of its vacuum expectation value, v v . This means we 
can choose the scale such that the Yukawa coupling h v 
is of the same order as all the other electroweak Yukawa 
coupling constants; in this case, if v v ~ 10 GeV and 
h v ~ 100, we get the correct scale for M v . 

Another example we give is one in which the Dirac 
neutrino masses are the same order as the top quark 
mass, and the Majorana mass scales are of the order 
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M GUT = 10 16 GeV, albeit with h lv mildly tuned. For 
this choice we have, 



VIII. RESOLVING THE VACUUM REALIGN- 
MENT PROBLEM VIA SPLITTING 



J177 



25.000, 



h p = 328.95, h' = 208.26, h" = 793.17, 



A p y=-i — ^4fl=3,y= 
B p y=-\ = -B_R=3,y= 



.1 = 900.00, 
-1 = 600.00, 



which gives 



If we then set 



Throughout this paper our analysis has depended on a 
particular choice of alignments for the vacuum expecta- 
(66) tion values of the two Aj-triplet Higgs fields localized to 
the wall. As per usual, finding valid VEVs for these fields 
involves finding global minima for the full scalar potential 
of the theory, which we have put in Appendix [A] For the 
gauge charged triplet, p, we assigned a VEV of the form 



m'p = 174 GeV. 



M = 8.65 x 10 15 GeV, M v = 6.85 x 10 15 GeV, (68) 

we get the masses of the left neutrino mass eigenstates 
for this parameter choice to be the same as those for 
the first parameter choice in Eqn. [65] Given that the 
M v = 6.85 x 10 15 GeV, if we have v v ~ 10 13 GeV, then 
we obtain h v ~ 100 as desired. 

In summary, we have shown that there exist param- 
eter choices within the model such that the fermion 
mass hierarchy, light neutrinos and the correct neutrino 
mass squared differences, quark mixing and a tribimax- 
imal lepton mixing matrix are reproduced from a set of 
Yukawa coupling constants which are of the same order 
of magnitude. In light of the recent results of the Daya 
Bay and RENO neutrino experiments |3] , which found 
that 7.9° < #13 < 9.6° for the PMNS matrix, exact tribi- 
maximal mixing is now excluded. However, a number of 
assumptions in this analysis were made which ensured an 
exact tribimaximal form, namely that m p = m'p = m!' 
and requiring that the left electron diagonalization ma- 
trix was close to trivial. It is obvious that deviations from 
tribimaximal mixing will occur if these assumptions are 
relaxed, and it is clear that the correct #13 mixing angle 
can also be generated in this model. 



The 



assumption generically results 



in a normal neutrino mass hierarchy, since the only 
possible hierarchies for this parameter choice are either 
m i < m 2 < "^3 j which is the normal hierarchy, or 
mi > 777.2 > TJI3 which is phenomenologically unaccept- 
able. Breaking this assumption can then lead to inverted 
or quasidegenerate neutrino mass hierarchies as well as 
deviations from tribimaximal mixing. 

Like many models based on A4, our analysis relied on 
the alignments of the vacuum expectation values of two 
^-triplet Higgs fields being different. However, when in- 
teractions between the two flavons are switched on, this 
arrangement is generally destroyed and the VEVs of the 
two fields align, leading to the vacuum realignment prob- 
lem. We will discuss the resolution of this problem within 
our model in the next section. 



(v 



), which induces a breaking of A4 — » Z3, and 



pi u pi v p 

for the neutrally charged triplet flavon p a VEV of the 
(67) form (0,v v ,0) which breaks A4 down to Zq. After these 
fields are localized and gain tachyonic masses on the wall, 
one can show that these VEV patterns are indeed valid 
global minima of the respective self-interaction potentials 

However, when one 



V p and V v given in Eqs. A2 and A3 



turns on interactions between the Higgs triplets and the 
1, V and 1" quintet scalars, this is no longer guaranteed 

[lain]- 

The problem is that once the cross-talking interactions 
are switched on, minimization of the potential yields a 
larger number of independent equations than known vac- 
uum expectation values (v p , v^, (4> w ), ( ( f ) ' w )i (^w))- This 
means that enforcing the desired vacuum alignment re- 
quires an unnatural fine-tuning of the parameters of the 
scalar potential. The troublesome terms are those which 
generate soft A4 — > Z2 breaking mass terms in the po- 
tential for p w after tp attains its VEV, and likewise gen- 
erate soft A4 — > Z3 breaking terms for the potential for 
tpo when p w condenses. These unwanted interactions in- 
volve coupling 1 , 1", 3 S and 3 Q products of p and p, for 
example (p^p)y (ipip)i" [22 . In addition to interactions 
with (p, there are other interactions generating A4 — > Z2 
breaking mass terms due to the presence of the 1' and 
1" Higgs quintets $' and <&", for example (p 1 p) V i&<&' . 
Analogous interactions involving just ip, $' and $" (and 
technically $ as well) are not problematic since ip, $' 
and <i>" all break A4 to the same Z2 subgroups, although 
there are interactions such as $"'(p<p)i" which provide 
corrections to the potential for ip as well as that for p. 

In our model, we follow the approach in 23 25 and 
split the triplet flavons in the extra dimension to ensure 
the alignment. Typically, if we suppress the troublesome 
interactions such that their effective interaction strengths 
are extremely small, we simply expect a small classical 

for 



\pv p 



correction to the desired vacua for p w of order 
quartic interactions of p and ip for instance, and sim- 
ilarly for other troublesome mass terms impacting the 
alignment of <p. Since p is a quintet under SU(5), the 
electroweak Higgs doublet p w arising from it is in gen- 
eral displaced from y = 0. Since the singlet flavon p is 
always localized at y — 0, it is in principle possible to sep- 
arate the flavons sufficiently so that the operators in V pip 
are naturally suppressed. As previously mentioned, we 
also have troublesome interactions involving p and the 1' 
and 1" Higgs quintets. Thus we place the localized mode 
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of p, p w , on the opposite side of the wall to those of $, 
<fr' and The interactions coupling all three types of 
non-trivially ^-charged flavons such as §"^(p<f)yi are 
naturally suppressed by the splitting of the A4 — » Z3 
breaking sector generated by p from the A± — > Z2 break- 
ing sector generated from cp, and 



4r 



I I 
I I 
I I 
I I 
I I 
I 1 
I I 
I I 
I I 
I I 
I I 
I I 
I 1 
I 1 
I 1 
I I 
I I 
I I 
f \ 



3- 



FIG. 7: Profiles for the localized triplet flavons p w and 
<po as well as those for the 1, 1' and 1" Higgs fields 4>w , 
<p' w and ^ for the parameter choices with A p y = -i and 



B p y=-i as given in Equation 62 the same localization 



parameters given for <£>, $' ana<£ > '' chosen in Sec. VII 
and d = 500 for if. 

In our example analysis, for the first choice with 
A p y=-i and B p y=~\ given in Eq. |62[ the same local- 
ization parameters for $, $' and $ 



VII 



chosen in Sec 

we choose d, as defined in Eq. [31] and which determines 
how localized <po is to the wall, to be set to d — 500.00. 
A graph of the profiles for p w , ipo and those of <fi w , <f>' w , 
and for this parameter choice is shown in Fig. [7} 

The potentially most troublesome interactions are 
the cubic interactions of the form (p' p) 3s a .tp and 

(|fV)i,i',i"* fl=1,1 f° r which the effective 4D cou- 
plings are proportional to the integrals J p {y)p Vo (jj) dy 

and J Ppw{y)p<p (y)Pw,w',w"(y) dy respectively. One has 
to be careful in calculating the magnitude of the cor- 
rections from these interactions since quartic and cubic 
scalar couplings have different mass dimension; in 4+1D 
a quartic coupling has mass dimension —1 while a cubic 
coupling has mass dimension | . For a generic cubic inter- 
action of the types in consideration, if the cubic coupling 
is a, and the quartic self-couplings for a generic A4 triplet 
are of order A, provided that the cubic interaction pro- 
vides just a small perturbation, then we anticipate that 
this perturbation is proportional to f . If we choose these 

numbers to be natural, that is a = aA^ w and A = AA^^ 
where Aow is the UV cut-off for our theory and a and 

A are dimcnsionless numbers of order 1, this means that 

3 

j ~ Ap W . An overlap integral O3 of the profiles of the 
scalars involved in a cubic interaction is k^ times the 



non-dimensionalized integral O3, while that for a quartic 
interaction self-interaction O4 is kOi ~ k. If we now con- 
sider (p^ p) 3s a .ip, after p w attains a VEV, this interaction 
provides a correction linear in ipo and so the correction 

3 _i . v 2 

to v v is then of order Af~, w k 2 O3-5-, and likewise when 

tpo condenses this interaction provides a A4 — > Z2 mass 

3 x ~ 

term for p w which provides a correction A^ w k~^03^, 

where O3 — J Pp w {y)p Vo {y) dy. We do not know what 
exactly the scales k and Ajjw are, but for the purpose 
of this analysis we will assume the worst case scenario 
as far as these corrections are concerned: k ~ 1 TeV and 
Adw ~ Mpi anc k — 10 19 GeV. For this first parame- 
ter choice, J p 2 pw {y)Pip {y) dy is of the order of 10~ 35 , 
which yields the correction to (ipo) from this interaction 
to be of order lO -7 ^, and the correction to (p w ) to be 
of order 10~ la v p . Likewise, for the interactions of the 
form (pVji.i'.i''^^ 1 ' 1 i the relevant overlap integral 
0-3 = 1 Ppw(y)p lpo (y)Pw,w',w"(y) dy is of order 10~ 31 and 
after the non-triplet gauge charged Higgs fields attain 
VEVs these interactions give corrections linear in ipo and 
p w , which turn out to be of the order of 10~ 3 v v and 
lO -8 ^ to (<po) and (p w ) respectively. 

With regard to quartic cross-talk interactions, we do 
not have to worry about scaling of the coupling constants 
for obvious reasons. For the quartic interactions coupling 
p and tp, the effective 4D quartic couplings are all pro- 
portional to the integral / p 2 {y)p^ (y) dy which for this 
first parameter choice turns out to be of order 10 -53 , 



yielding corrections of order (10 ^)v i 



10 



-o5 



v p to 



( Pw ) and (lO" 53 ^)^ 



10 



-52, 



to (ip ) 



Next, we deal with the quartic interactions cou- 
pling p, $ and one of $' or which are of the 
form (p' p)r=i'.i"^^ r=1 ,x . After dimensional re- 
duction the effective 4D coupling constants for these 
interactions are proportional to the overlap integral 
/ Ppw(y)Pw(y)Pw',w"(y) dy. For this parameter choice 
these overlap integrals are of order 10 -22 , yielding cor- 
rections of order (10~ 22 M^K ),, _ 10- 22 w p to 

P 

(Pw)- 

Finally, there are quartic interactions coupling both 
p and if to one of <&, $' or which are of the form 
[(<fip) 3s .p^] R $ R *. To leading order, since ((<po) (<Po)hs = 



((0,^,0).(0,^,0)) 3s 



0, this interaction does not 



give a leading order correction to (p w ), although there 
are corrections for (fo) . After dimensional reduction, 
the effective 4D coupling strength is proportional to 
the overlap integral f p^ (y)p pw (y)Pw,w',w"(y) dy, and 
for the first parameter choice this overlap integral is 
of the order 10 -39 , leading to corrections of order 



(10 



-39 v p{<f>. 



i?=i,i / ,i / ' 



)v ip ~ 10 3 v v to ((fa). Obviously any 



Jin v<p 

second order corrections to (p w ) will be much smaller 
than this. 

For the second choice with parameters given in Eq. [66] 
with d = 900.00, we get the value of the relevant over- 
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lap integrals to be of order 10 50 for J p 2 pw {y)p Vo {y) dy, 

lCr 42 for J Pp W (y)p ipo (y)p w , W ', W "(y) dy, icr 77 for 
f p 2 pw(y)pl (y) d v> 10 ~ 19 for I Ppw(y)p 2 u(y)pw,w"(y) dy, 

and 1(T 54 for / Pp W (y)pfy (y)p w ,w',w''(y) dy. All the cor- 
rections to (tpa) are not only suppressed by the over- 
lap integrals but by either ^f- ~ 10~ 12 or its square 

(^) 2 ~ 10~ 24 , hence doing a similar analysis to that 
above shows that the corrections to (ipo) are extremely 
negligible. Doing a similar analysis as above shows that 
the the overlap integrals are small enough to overcome 
the ratio ~ ~ 10 12 or its square (^) 2 ~ 10 24 (whichever 

P 3 _i 

is relevant) as well as the ratio Ap W k 2 in the case of 

the cubic interactions and still ensure that the correc- 
tions to (p w ) are more than several orders of magnitude 
less than v p . 

IX. CONCLUSION 

We found that the A A extension of the SU(5) 4+1D 
domain- wall braneworld model of |16j can generate large 
mixing angles in the lepton sector. We explicitly demon- 
strated parameter values that yield tribimaximal lepton 
mixing with a normal neutrino mass hierarchy together 
with successful predictions for the hierarchical charged 
fermion masses and quark mixing angles. Through small 
departures from this parameter point, the small but 
nonzero #13 leptonic mixing angle can be generated, and 
the neutrino mass spectrum altered to give an inverted 
hierarchy or a quasi-degenerate pattern. This is a signif- 
icant extension of the results found in [17] . 

We also discovered that the troublesome interactions 
which are responsible for the vacuum realignment prob- 
lem in analogous 4D models could be suppressed by split- 
ting the profile of p, which initiates A4 — > Z3 breaking, 
from the profiles of the A4 — > I2 breaking flavons (p, 
and This led to exponential suppression of the 
overlap integrals of these profiles determining the effec- 
tive 4D coupling constants for interactions between the 
A-triplets as well as the troublesome interactions which 
involve <&' and ensuring that the contributions of 
these interactions to the vacua of the localized compo- 
nents of p and tp were sufficiently small compared to the 
classical vacua of their respective self-interaction poten- 
tials. This maintained the desired vacuum alignments 
for these flavons required to generate large lepton mixing 
angles. 

We have shown in this paper and in |17j that domain- 
wall braneworld models with an SU (5) gauge group are 
generically good for generating the fermion mass hierar- 
chy and quark mixing, and that with the inclusion of the 
discrete flavor group A4 we can attain large lepton mix- 
ing angles. Further work along these lines could involve 
extending either the gauge group (to 50(10) for example 
[37] . or to even larger groups such as E 6 [38]), the dis- 
crete flavor group (for example, to T' which is the double 
cover of At), or both. Further work could also look at the 



case where gravity is included or quantum corrections to 
the fermion mass spectra. 



Acknowledgments 

This work was supported, in part, by the Australian 
Research Council and the Commonwealth of Australia. 
We would also like to thank Ferruccio Feruglio for use- 
ful discussions as well as the Dipartimento di Fisica 'G. 
Galilei' at the Universita di Padova for its warm hospi- 
tality. 



Appendix A: The Higgs Flavon Scalar Interaction 
Potential 

In this appendix, we give the potentials describing the 
interactions amongst different Higgs flavon fields local- 
ized to the domain wall. The background scalar field 
potential V vx yielding the background kink-lump solu- 
tion was given in Eqn. [2] in Sec. [Tl] The potentials cou- 
pling the flavons to the background fields rj and x, W$, 
Wq,it, W p and W v , were given in Sec. [Vjin Eqns. 
[20] and [29j which after localization of the scalars and 
an appropriate choice of parameters contribute tachy- 
onic masses to the effective self-interaction potentials 
for the lightest localized modes of these fields. Thus 
let Vp, V v , V$, V$/, and V$" be potentials containing 
the quartic self-interactions for each of these fields along 
with their localization potentials. Let the potentials V pv> , 
Vp$$'$", V v $$'$", V pv ^w, and V$$<$» be those con- 
taining cross-talk interactions between the fields in the 
subscripts. Then the full scalar potential of the theory, 
V, is given by 

V = V nx + V P + V V + V* + V*. + V*« + V pv 
+ Vp<j>$'$" + V v $$'$" + V P cp$w + V$$'$", 

where 

V P = A p (p f p)i(pV)i + \ 2 p {p ] p)v{J p)v 

+ \ p [( P V)3,.(p t p) 3s ] l + Pha.{p ] P)*a\ ! 

+ iX P [(p i phs.(p i pha] 1 + Wp 

= A>V) 1 (pV)i + A>V)i'(P t P)i» 

+ A P [(PV) 3 ,VP) 3S ] l + A 4 [( P V) 3 aVp) 3 a] x 

+ *Aj [(p t p) 3 s-(p t pha\ l + PpVp)i + Kn^PW 
+ 2\ pxl (p^p) 1 Tr( X 2 ) + \ PX 2{P ] (x T Yp)i 
+ V)x(/ ,t X T p)i»7 ; 

(A2) 
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V, 



+ A 2 (w)i'(w)i" + ^ [(w)3 S (w)3 S ]! +W V 
= ^(<PW)i + ^(w)i(w)i 

+ A£(w)l'(v?V)l" + ^[(^3^(^)38]! 

+ *4(w)i + Kn(PV)iV 2 + 2A ¥ , x (^)iTr(x 2 ), 

(A3) 



= A $ «(($*)t<^) 2 + 4 s ($«)t$ fl + A^„($«)t$V 
+ 2A^ xl ($«)tcf>«rr( X 2 ) + A $ , x2 ($«)t( x r ) 2 cf>« 
+ \*R nx {* R ?X T * R ri, for J2 = 1, 1', 1", 

(A4) 



= Ai $ ,($t$)($'t$') + A| v ($t$')($t$) 
+ A^„($ t $)($" t $") + A| 4 „($t$")($"t$) 
+ A|„i,„($' t $ / )($" t $") + A|„ I) „($' t $")($" t $') 
+ Ai»'*»(* t * , )(* t *") + AiW'^X*"'*) 
+ A| $ , <E> „($ / t$)($'t$") + AlV*"^*')^*') 
+ A| $ , <E> „($"t$)($"t$') + AlW'^'X*' 1 *") 

(A9) 



+ A p V (p f P)i(w)i + ^(pMi'Cw)!" 

+ A 2 ;(p f p)i"(w)i' + aI. \(p^p) 3s .( w ) 3s 



+ iX 



(A5) 



W*" = A^(pt p )!($t$) + A p */(pt p ) 1 ($t#') 

+ A p * < Vp) 1 ($t$') + A p **»(ptp) 1 „($t$') 

+ a; m „(pV)i» + A^VpV^") 

+ AjJ < *»(pV)i'(*" t * , ) + A^[(p t *).(p t *)] 1 
+ \% [($tp).($t p) ] x + A 2^„ [(^.(ptd,")] i 

+ A 2 ^ 4) „[(<i>'tp).(cf>"t p )] i 

(A6) 

= A v $($ t $)(( / 9( / 5)l + A^' ($'*$') (W>)l 

+ A ¥ ,«»($"t$")(w)i + V*'( $t *')(w)i« 
+ A; $$ ,($ /t $)( W )^ + A v **» (* t $")(w)i' 
+ A; $$ „($"t$)( w ) 1 „ + A^^($'t$'')( w ) 1 „ 

+ A; 4 ,*„(*"t$ , )(w)i'. 

(A7) 

+ <W'(w f )i'*" + 5* P ^»&'Hptp)i" 

+ [( w ) 3s .pt] ^ + A; v$ $ f [p.(<p<p) 3 J x 

+ \ PV *> [( W ) 3s .p f ] + A^*,** [p.(w) 3 J r 

+ a pv $» [(w) 3s .p f ] + a;^,,*" 1 [p-(w) 3s ] lf „ 

(A8) 
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